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Abstract. PGA (ProGram Algebra) is an algebra of programs which 
concerns programs in their simplest form: sequences of instructions. 
Molecular dynamics is a simple model of computation developed in the 
setting of PGA, which bears on the use of dynamic data structures in pro- 
gramming. We consider the programming of an interpreter for a program 
notation that is close to existing assembly languages using PGA with the 
primitives of molecular dynamics as basic instructions. It happens that, 
although primarily meant for explaining programming language features 
relating to the use of dynamic data structures, the collection of primi- 
tives of molecular dynamics in itself is suited to our programming wants. 
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1 Introduction 

In this paper, we consider the programming of an interpreter for a program nota- 
tion that is close to existing assembly languages using PGA (ProGram Algebra). 
With that we carry on the line of research with which a start was made in [3]. 
The object pursued with that line of research is the development of a theoret- 
ical understanding of possible forms of sequential programs, starting from the 
simplest form of sequential programs, and associated ways of programming. The 
view is taken that sequential programs in the simplest form are sequences of 
instructions. PGA, an algebra of programs in which programs are looked upon 
as sequences of instructions, is taken for the basis of the development aimed at. 
Until now, the research was rather focussed on the development of a theoretical 
understanding of possible forms of sequential programs. The work presented in 
the current paper is primarily concerned with programming using the simplest 
form of sequential programs and molecular dynamics. 



* This research was partly carried out in the framework of the Jacquard-project Sym- 
biosis, which is funded by the Netherlands Organisation for Scientific Research 
(NWO). 



For the programming of the interpreter, we use PGA with the primitives of 
molecular dynamics as basic instructions. Molecular dynamics is a simple model 
of computation bearing on the use of dynamic data structures in programming. 
In this model, states of computations resemble collections of molecules composed 
of atoms and computations take place by means of actions which transform the 
structure of molecules like in chemical reactions. Molecular dynamics has been 
developed mainly in the setting of PGA. The model introduced in the current 
paper is a small expansion of the model that was first described informally in [1] 
and formally in [5]. 

In the line of research carried on, the view is taken that the behaviours ex- 
hibited by sequential programs on execution are threads as considered in BTA 
(Basic Thread Algebra).^ A thread proceeds by doing steps in a sequential fash- 
ion. A thread may do certain steps for having itself affected by some service 
or affecting some service;. In the current paper, we will use a generalization of 
the use mechanism introduced in [8] and a complementary mechanism of that 
mechanism for such interactions between threads that are behaviours exhibited 
by sequential programs on execution and services that deal with steps that relate 
to molecular dynamics. A slightly different generalization of the use mechanism 
from [8] is introduced in [6] under the name action transforming thread- service 
composition. 

A hierarchy of program notations rooted in PGA is introduced in [3]. In the 
current paper, we consider the programming of an interpreter for one program 
notation that belongs to this hierarchy. The program notation in question, called 
PGLD (ProGramming Language D), is a very simple program notation which 
is close to existing assembly languages. The hierarchy also includes a program 
notation, called PGLS (ProGramming Language for Structured programming), 
that supports structured programming by offering a rendering of conditional 
and loop constructs instead of (unstructured) jump instructions. Each PGLS 
program can be translated into a semantically equivalent PGLD program. In [7], 
a variant of PGLD with indirect jump instructions is introduced. In the current 
paper, we show how to adapt the interpreter for PGLD to the presence of indirect 
jump instructions. 

This paper is organized as follows. First, we review BTA, PGA, and PGLD 
(Sections 2, 3, and 4). Next, we extend BTA with the mechanisms for interaction 
between threads and services used, introduce a state-based approach to describe 
services, and give a state-based description of services for molecular dynamics 
(Sections 5, 6, and 7). Following this, we give PGA programs for creating repre- 
sentations of PGLD programs by molecules and a PGA program for interpreting 
those representations (Sections 8 and 9). After that, we introduce the variant 
of PGLD with indirect jump instructions and adapt the above-mentioned PGA 
programs to the presence of indirect jump instructions (Sections 10 and 11). 
Finally, we make some concluding remarks (Section 12). 

^ In [3], basic thread algebra is introduced under the name basic polarized process 
algebra. Prompted by the development of thread algebra [8] , which is a design on 
top of it, basic polarized process algebra has been renamed to basic thread algebra. 
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Table 1. Axiom of BTA 



a; < tau > J/ = a; < tau > a; Tl 



2 Basic Thread Algebra 

In this section, we review BTA, a form of process algebra which is tailored 
to the description of the behaviour of deterministic sequential programs under 
execution. The behaviours concerned are called threads. 

In BTA, it is assumed that there is a fixed but arbitrary finite set of basic 
actions A with tau ^ A. We write ^tau for A U {tau}. The members of ^tau are 
referred to as actions. 

The intuition is that each basic action performed by a thread is taken as a 
command to be processed by a service provided by the execution environment of 
the thread. The processing of a command may involve a change of state of the 
service concerned. At completion of the processing of the command, the service 
produces a reply value. This reply is either T or F and is returned to the thread 
concerned. 

Although BTA is one-sorted, we make this sort explicit. The reason for this 
is that we will extend BTA with an additional sort in Section 5. 

The algebraic theory BTA has one sort: the sort T of threads. To build terms 
of sort T, BTA has the following constants and operators: 

— the deadlock constant D : T; 

— the termination constant S : T; 

— for each a S Aau, the binary postconditional composition operator _ < a > _ : 
T X T ^ T. 

Terms of sort T arc built as usual (see e.g. [13, 14]). Throughout the paper, we 
assume that there are infinitely many variables of sort T, including x, y, z. 

We use infix notation for postconditional composition. We introduce action 
prefixing as an abbreviation: a o p, where p is a term of sort T, abbreviates 
p<a>p. 

Let p and q be closed terms of sort T and a G Aau- Then p <a>q will 
perform action a, and after that proceed as p if the processing of a leads to the 
reply T (called a positive reply), and proceed as q if the processing of a leads 
to the reply F (called a negative reply). The action tau plays a special role. It 
is a concrete internal action: performing tau will never lead to a state change 
and always lead to a positive reply, but notwithstanding all that its presence 
matters. 

BTA has only one axiom. This axiom is given in Table 1. Using the abbrevia- 
tion introduced above, axiom Tl can be written as follows: x < tau >y = tauox. 

Each closed BTA term of sort T denotes a finite thread, i.e. a thread of which 
the length of the sequences of actions that it can perform is bounded. Guarded 
recursive specifications give rise to infinite threads. 
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Table 2. Axioms for guarded recursion 



{X\E) = {tx\E) if X = tx e E RDP 
E ^ X = {X\E) if X € V{E) RSP 



A guarded recursive specification over BTA is a set of recursion equations 
E = {X ^ tx \ X E V}, where F is a set of variables of sort T and each tx 
is a term of the form D, S or t <a>t' with t and t' BTA terms of sort T that 
contain only variables from V. We write V(S) for the set of all variables that 
occur on the left-hand side of an equation in E. We are only interested in models 
of BTA in which guarded recursive specifications have unique solutions, such as 
the projective limit model of BTA presented in [2] . A thread that is the solution 
of a finite guarded recursive specification over BTA is called a finite-state thread. 

We extend BTA with guarded recursion by adding constants for solutions 
of guarded recursive specifications and axioms concerning these additional con- 
stants. For each guarded recursive specification E and each X G V(i?), we add a 
constant of sort T standing for the unique solution of E for X to the constants 
of BTA. The constant standing for the unique solution of E for X is denoted by 
{X\E). Moreover, we add the axioms for guarded recursion given in Table 2 to 
BTA, where we write {tx\E) for tx with, for all Y G V(£'), all occurrences of 
Y in tx replaced by {Y\E). In this table, X, tx and E stand for an arbitrary 
variable of sort T, an arbitrary BTA term of sort T and an arbitrary guarded re- 
cursive specification over BTA, respectively. Side conditions are added to restrict 
the variables, terms and guarded recursive specifications for which X, tx and E 
stand. The equations {X\E) = {tx\E) for a fixed E express that the constants 
{X\E) make up a solution of E. The conditional equations E ^ X = {X\E) 
express that this solution is the only one. 

We will use the following abbreviation: a'^, where a G A^u, abbreviates 
{X\{X = aoX}). 

We will write BTA+REC for BTA extended with the constants for solutions 
of guarded recursive specifications and axioms RDP and RSP. 

In [4], we show that the threads considered in BTA-I-REC can be viewed as 
processes that are definable over ACP [10] . 

Closed terms of sort T from the language of BTA+REC that denote the 
same infinite thread cannot always be proved equal by means of the axioms 
of BTA+REC. We introduce the approximation induction principle to remedy 
this. The approximation induction principle, AIP in short, is based on the view 
that two threads are identical if their approximations up to any finite depth are 
identical. The approximation up to depth n of a thread is obtained by cutting 
it off after performing a sequence of actions of length n. 

AIP is the infinitary conditional equation given in Table 3. Here, following [3], 
approximation of depth n is phrased in terms of a unary projection operator 
7r„ : T ^ T. The axioms for the projection operators are given in Table 4. In 
this table, a stands for an arbitrary member of ^tau- 
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Table 3. Approximation induction principle 



t\n>a^r,{x) = 7r„(j/) ^ x = y AIP 

Table 4. Axioms for projection operators 

7ro(a;) = D PO 

7r„+i(S) = S PI 

7r„+i(D) = D P2 

7r„+i(a; < a> J/) = 7r„(a;) <o> 7r„(y) P3 



We will write BTA+REC+AIP for BTA+REC extended with the projection 
operators and the axioms from Tables 3 and 4. 



3 Program Algebra 

In this section, we review PGA, an algebra of sequential programs based on 
the idea that sequential programs are in essence sequences of instructions. PGA 
provides a program notation for finite-state threads. 

In PGA, it is assumed that there is a fixed but arbitrary finite set 21 of basic 
instructions. PGA has the following primitive instructions: 

— for each a G 21, a plain basic instruction a; 

— for each a S 2t, a positive test instruction +a; 

— for each a e 21, a negative test instruction —a; 

— for each Z G N, a forward jump instruction ij^l; 

— a termination instruction !. 

We write 3 for the set of all primitive instructions. 

The intuition is that the execution of a basic instruction a may modify a 
state and produces T or F at its completion. In the case of a positive test in- 
struction +a, basic instruction a is executed and execution proceeds with the 
next primitive instruction if T is produced and otherwise the next primitive 
instruction is skipped and execution proceeds with the primitive instruction fol- 
lowing the skipped one. In the case where T is produced and there is not at least 
one subsequent primitive instruction and in the case where F is produced and 
there arc not at least two subsequent primitive instructions, deadlock occurs. 
In the case of a negative test instruction —a, the role of the value produced is 
reversed. In the case of a plain basic instruction a, the value produced is disre- 
garded: execution always proceeds as if T is produced. The effect of a forward 
jump instruction is that execution proceeds with the Z-th next instruction of 
the program concerned. If I equals or the /-th next instruction does not exist, 
then #Z results in deadlock. The effect of the termination instruction ! is that 
execution terminates. 

PGA has the following constants and operators: 
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Table 5. Axioms of PGA 



{x\y)\z = x;{y;z) PGAl 
{x'^Y = x'^ PGA2 
a;" ; y = a;" PGA3 
{x;yY = x;{y\ a:)" PGA4 



— for each u G J, an instruction constant u ; 

— the binary concatenation operator _ ; _ ; 

— the unary repetition operator . 

Terms are built as usual. Throughout the paper, we assume that there are in- 
finitely many variables, including x,y, z. 

We use infix notation for concatenation and postfix notation for repetition. 

Closed PGA terms are considered to denote programs. The intuition is that 

a program is in essence a non-empty, finite or periodic infinite sequence of prim- 
itive instructions.^ These sequences are called single pass instruction sequences 
because PGA has been designed to enable single pass execution of instruction 
sequences: each instruction can be dropped after it has been executed. Programs 
are considered to be equal if they represent the same single pass instruction se- 
quence. The axioms for instruction sequence equivalence are given in Table 5. 
In this table, n stands for an arbitrary natural number greater than 0. For 
each n > 0, the term x" is defined by induction on n as follows: = x and 
a;"+i = x-jX"". The unfolding equation x'^ = a;;x" is derivable. Each closed PGA 
term is derivably equal to a term in canonical form, i.e. a term of the form P or 
P; Q'^ , where P and Q are closed PGA terms that do not contain the repetition 
operator. 

Each closed PGA term is considered to denote a program of which the be- 
haviour is a finite-state thread, taking the set 21 of basic instructions for the set A 
of actions. The thread extraction operator |_| assigns a thread to each program. 
The thread extraction operator is defined by the equations given in Table 6 (for 
a G 2t, Z e N and u <E 3) and the rule given in Table 7. This rule is expressed 
in terms of the structural congruence predicate _ — _, which is defined by the 
formulas given in Table 8 (for n, to, ? G N and Wi, . . . , m„, wi, . . . , v„i+i G 3). 

The equations given in Table 6 do not cover the case where there is a cyclic 
chain of forward jmups. Programs are structural congruent if they are the same 
after removing all chains of forward jumps in favour of single jumps. Because 
a cyclic chain of forward jumps corresponds to ^0, the rule from Table 7 can 
be read as follows: if x starts with a cyclic chain of forward jumps, then \x\ 
equals D. It is easy to see that the thread extraction operator assigns the same 
thread to structurally congruent programs. Therefore, the rule from Table 7 can 
be replaced by the following generalization: x = y => |a;| = |y|. 

A periodic infinite sequence is an infinite sequence with only finitely many subse- 
quences. 
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Table 6. Defining equations for thread extraction operator 



a\ = aoD \#l\^D 

a;x\ = ao\x\ |#0 ; a;| = D 

+a|=aoD |#l;a-l = |xl 

+a;x| = ^a> |#2;a;| |#; + 2 ; ti| = D 

-a\ = aoD \#l + 2;u;x\ = \#l + l;x\ 

-o ;a;| = |#2 ;a;| <a> |a;| |!| = S 

|!;a:|=S 

Table 7. Rule for cyclic jump chains 



x^#0;y ^ \x\ = D 

Table 8. Defining formulas for structural congruence predicate 

#n + 1 ; wi ; . . . ; w„ ; #0 = #0 ; Ml ; . . . ; w„ ; #0 

#n + 1 ; Ml ; . . . ; M„ ; #m = #m + n + 1 ; mi ; . . . ; w„ ; #m 

(#n + / + 1 ; Ml ; . . . ; = (#Z ; Mi ; . . . ; w„)" 

#m + n + Z + 2 ; Ml ; . . . ; M„ ; (vi ; . . . ; Vm+i)'^ = 

#n + Z + 1 ; Ml ; . . . ; M„ ; (mi ; . . . ; Mra+i)"^ 

X = 2: 

xi = yi Ax2 = y2 =^ xi ; X2 = yi ■,y2 A xi^ = yx" 



Let i5 be a finite guarded recursive specification over BTA, and let Px be a 
closed PGA term for each X G V(i?). Let E' be the set of equations that results 
from replacing in E all occurrences of X by \Px\ for each X e V(£'). If E' can 
be obtained by applications of axioms PGA1-PGA4, the defining equations for 
the thread extraction operator and the rule for cyclic jiunp chains, then \Px\ is 
the solution of E for X. Such a finite guarded recursive specification can always 
be found. Thus, the behaviour of each closed PGA term, is a thread that is 
definable by a finite guarded recursive specification over BTA. Moreover, each 
finite guarded recursive specification over BTA can be translated to a closed 
PGA term of which the behaviour is the solution of the finite guarded recursive 
specification concerned (cf. Section 4 of [12]). 

Closed PGA terms are loosely called PGA programs. PGA programs in which 
the repetition operator do not occur are called finite PGA programs. 

4 The Program Notation PGLD 

In this section, we review a program notation which is rooted in PGA. This 
program notation, called PGLD, belongs to a hierarchy of program notations 
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introduced in [3]. PGLD is close to existing assembly languages. It has absolute 
jump instructions and no explicit termination instruction. 

In PGLD, like in PGA, it is assumed that there is a fixed but arbitrary finite 
set of basic instructions 21. Again, the intuition is that the execution of a basic 
instruction a may modify a state and produces T or F at its completion. 

PGLD has the following primitive instructions: 

— for each a £ 21, a plain basic instruction a: 

— for each a G 21, a positive test instruction +a; 

— for each a e 21, a negative test instruction —a; 

— for each Z G N, a direct absolute jump instruction 

PGLD programs have the form u\;...;Uk, where Ui,...,Uk are primitive in- 
structions of PGLD. 

The plain basic instructions, the positive test instructions, and the negative 
test instructions arc as in PGA. The effect of a direct absolute jump instruction 
##i is that execution proceeds with the /-th instruction of the program con- 
cerned. If ##Z is itself the Z-th instruction, then deadlock occurs. If I equals 
or I is greater than the length of the program, then termination occurs. 

We define the meaning of PGLD programs by means of a function pglci2pga 
from the set of all PGLD programs to the set of all PGA programs. This function 
is defined by 

pgld2pga(ui ; . . . ; Ufe) = (V'i(wi) ; • • • ; V'fe(^ifc) ; ! ; 0'^ > 

where the auxiliary functions ij^j from the set of all primitive instructions of 
PGLD to the set of all primitive instructions of PGA are defined as follows 
(1 < j <fc): 

= if.?<;<fc, 

^j{##l) = #k + 2-{j-l) \fO<l<j, 

V)^.(##Z) = ! ifZ = OVZ>fc, 

V'j (^) = ^ if u is not a jump instruction . 

Let P be a PGLD program. Then pgld2pga(P) represents the meaning 
of P as a PGA program. The intended behaviour of P is the behaviour of 
pgld2pga(P). That is, the behaviour of P, written |P|pgld, is |pgld2pga(P)|. 

We use the phrase projection semantics to refer to the approach to semantics 
followed in this section. The meaning function pgld2pga is called a projection. 

5 Threads Interacting with Action Transforming Services 

A thread may perform certain actions for having itself affected by some service 
or affecting some service. When processing an action performed by a thread, 
an action transforming service affects that thread in one of the following ways: 
(i) by returning a reply value to the thread at completion of the processing 
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of the action performed by the thread; (ii) by turning the processed action into 
another action. In this section, we introduce a mechanism that allows for services 
to affect threads in either way. The mechanism is a generalization of the use 
mechanism introduced in [8].^ We also introduce a complementary mechanism 
of that generalized use mechanism and an abstraction mechanism. The difference 
between the generalized use mechanism and the complementary mechanism is 
a matter of perspective: the former is concerned with the effect of services on 
threads and therefore produces threads, whereas the latter is concerned with 
the effect of threads on services and therefore produces services. The abstraction 
mechanism serves for concealment of the presence of internal actions, which arise 
primarily from the generalized use mechanism. 

We will use the generalized use mechanism and the complementary mech- 
anism of that mechanism for interactions between threads that are behaviours 
exhibited by sequential programs on execution and services that process actions 
that relate to molecular dynamics. 

It is assumed that there is a fixed but arbitrary finite set of foci T and a fixed 
but arbitrary finite set of methods M.. Each focus plays the role of a name of a 
service provided by the execution environment that can be requested to process 
a command. Each method plays the role of a command proper. For the set A 
of actions, we take the sot {/.m | / G .F, m S M}. Performing an action f.m is 
taken as making a request to the service named / to process command m. 

We introduce yet another sort: the sort S of services. S is considered to 
stand for the set of all services. We identify services with pairs {Hi,H2), where 
Hi:M~^ {T, F, M, B} and H2 :A4~^ — > Atau, satisfying the following conditions: 

Vm e M • 

(3a G M* • Hi{a ^ (m)) = M ^ Va' € M* ' Hi{a' ^ (m)) ^ {T, F}) , 
Va G M+,m G M • {Hi{a) = B Hi{a ^ (m)) = B) , 
Va G • {Hi{a) 7^ M <^ H2{a) = tau) . 

Let iJ be a service, and let Hi and H2 be the unique functions such that H = 
{Hi, H2). Then wc write rf{H) and af{H) for Hi and H2, respectively. Given a 
service H and a method m G ^A, the derived service of H after processing m, 
written £^H, is defined by rf{£^H){a) = rf{H){{m) -a) and af{£^H){a) = 
a/(F)((m}-a). 

A service H can be understood as follows: 

— if rf{H){{m)) = T, then the request to process m is accepted by the service, 
a positive reply is produced, m is turned into tau, and the service proceeds 

as -^H; 

am ' 

— if rf{H){{m)) = F, then the request to process m is accepted by the service, 
a negative reply is produced, m is turned into tau, and the service proceeds 



In later papers, the original use mechanism is also called thread-service composition. 
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— if rf{H){{m)) = M, then the request to process m is accepted by the service, 
no reply is produced, m is turned into af{H){{m)), and the service proceeds 

— if rf{H){{m)) = B, then the request to process m is rejected by the service. 
The three conditions imposed on services can be paraphrased as follows: 

— if it is possible that no reply is produced at completion of the processing of a 
command, then it is impossible that a positive or negative reply is produced 
at completion of the processing of that command; 

— after a request to process a command has been rejected, any request to 
process a command will be rejected; 

— a reply is produced at completion of the processing of a command if and 
only if the command is turned into tau. 

Wc introduce the following additional constant and operators: 

— the divergent service constant D : S; 

— for each f € J^, the binary use operator _//_:TxS— >T; 

— for each f € J^, the binary apply operator _«/_:TxS— >S. 

We use infix notation for use and apply. 

D is the unique service H with the property that rf{H){a) = B for all 
a & The operators - // - and - •/ - are complementary. Intuitively, p / f H 
is the thread that results from processing all actions performed by thread p 
that are of the form f.m by the service H. When an action of the form f.m 
performed by thread p is processed by the service H, that action is turned into 
another action and, if this action is tau, postconditional composition is removed 
in favour of action prefixing on the basis of the reply value produced. Intuitively, 
p»f H is the service that results from processing all basic actions performed by 
thread P that arc of the form f.m by service H . When an action of the form 
f.m performed by thread p is processed by service H, that service is changed 
into -j^H. 

am 

The axioms for the use and apply operators arc given in Tables 9 and 10. 
In these tables, / and g stand for an arbitrary foci from J^, m stands for an 
arbitrary method from M, and H is a variable of sort S. Axioms TSU3 and 
TSU4 express that the action tau and actions of the form g.m, where / =/= g, 
are not processed. Axioms TSU5-TSU7 express that a thread is affected by a 
service as described above when an action of the form f.m performed by the 
thread is processed by the service. Axiom TSU8 expresses that deadlock takes 
place when an action to be processed is not accepted. Axiom TSU9 expresses 
that the divergent service does not accept any action. Axiom TSA3 expresses 
that a service is not affected by a thread when the action tau is performed by 
the thread and axiom TSA4 expresses that a service is turned into the divergent 
service when an action of the form g.m, where f ^ g,is performed by the thread. 
Axioms TSA5 -TSA7 express that a service is affected by a thread as described 
above when an action of the form f.m performed by the thread is processed by 
the service. Axiom TSA8 expresses that a service is turned into the divergent 
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Table 9. Axioms for use operators 



S/fH = S 








TSUI 


D/fH = D 








TSU2 


tau o X / f H — tau o 


{x/fH) 






TSU3 


{x<g.m\>y) /f H 


= (x /fH)<g.mt>{y/fH) 


if 




TSU4 


{x<f.m>y)/fH 


= tauo(x/,^/f) 


if rf{H){{m)) = 


T 


TSU5 


{x<f.m>y)/fH 


= tBUo{y/f£^H) 


if rf{H){{m)) = 


F 


TSU6 


ix<f.m>y) /fH 


= {x<af{H){{m))>y)/f£^H 


if r/(//)((m)) = 


M 


TSU7 


(x<f.m\> y) If H 


= D 


if rf{H){{m)) = 


B 


TSU8 


{x < f.m > y) If D = 


= D 






TSU9 



Table 10. Axioms for apply operators 



S»f H = H 










TSAl 


D •/ -ff = D 










TSA2 


(tau o x) •/ H = X * 


fH 








TSA3 


[x < g.m >y)»fH 


= D 




if 




TSA4 


{x < f.m >y)»fH 


= X 4-H 

J am 




if r/(ff)((m)) 


= T 


TSA5 


{x < f.m \>y)»jH 


y *^ dm ^ 




if r/(if)((m)) 


= F 


TSA6 


{x < f.m>y) •f H 


= {x<af{H){{ 


m))>y) .f^H 


\frf{H)i{m)) 


= M 


TSA7 


{x < f.m> y) •f H 


= D 




if rf{H){(m)) 


= B 


TSA8 


(x < f.m > J/) •/ D 


= D 








TSA9 


iAn>0^ri{x)»f H = 


= D) ^ x*f H 


= D 






TSAIO 



service when an action performed by the thread is not accepted. Axiom TSA9 
expresses that the divergent service is not affected by a thread when an action 
of the form f.m is performed by the thread. 

We say that p Uf H is a divergent thread application if 7r„(p) •f H = D for 
ah n G N. Axiom TSAlO can be read as foUows: if p •/ is a divergent thread 
apphcation, then it equals D. 

The use operators introduced in [8] deals in essence with services H where 
af{H){a) = tau for all a S . For these services, the use operators introduced 
here coincide with those use operators. 

Let T stand for either BTA, BTA+REC or BTA+REC+AIP. Then we will 
write T + TSI for T, taking the set {f.m \ f G J^,m € M} for A, extended with 
the divergent service constant, the use and apply operators, and the axioms from 
Tables 9 and 10, with the exception of axiom TSAIO in the case where T does 
not stand for BTA+REC+AIP. 

The action tau is an internal action whose presence matters. To conceal its 
presence in the case where it does not matter after all, we also introduce the 
unary abstraction operator Ttau : T ^ T. 

The axioms for the abstraction operator are given in Table 11. In this table, 
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Table 11. Axiom for abstraction 



Ttau(S) = S TTl 

Ttau(D) = D TT2 

n,,{x<tau>y) =Tt,u{x) TT3 

rtau(a; < a > y) = Tuaix) <a> rtau(2/) TT4 



a stands for an arbitrary basic action from A. 

Abstraction can for instance be appropriate in the case where tau arises from 
turning actions of an auxiliary nature into tau on use of a service. An example 
of this case will occur in Section 9. Unlike the use mechanisms introduced in [8] 
and in the current paper, the use mechanism introduced in [9] incorporates 
abstraction. 

Let T stand for either BTA, BTA+REC, BTA+REC+AIP, BTA+TSI, 
BTA+REC+TSI or BTA+REC+AIP+TSI. Then we will write T+ABSTR for 

T extended with the abstraction operator and the axioms from Table 11. 

The equation rtau(tau") = D is derivable from the axioms of BTA+REC+ 
AIP+ABSTR. 

To simplify matters, from now on the set {f.m \ f £ J-,m £ A4} is taken as 
the set 21 of basic instructions when PGA or PGLD is concerned. Thereby no 
real restriction is imposed on the set 21: in the case where the cardinality of 
equals 1, all basic instructions have the same focus and the set A4 of methods 
can be looked upon as the set 21 of basic instructions. We use the convention 
to omit foci from PGA programs in which all basic instructions have the same 
focus. 

Strictly speaking, the propositions and theorems presented in this paper 
are proved in the algebraic theory obtained by: (i) combining PGA with 
BTA+REC+AIP+TSI+ABSTR, resulting in a theory with three sorts: a sort 
P of programs, a sort T of threads, and a sort S of services; (ii) extending the 
result by taking |_| for an additional operator from sort P to sort T and tak- 
ing the semantic equations and rule defining thread extraction for additional 
axioms. 

6 State-Based Description of Services 

In this section, we introduce the state-based approach to describe families of 
services that will be used in Section 7. This approach is similar to the approach 
to describe state machines introduced in [9]. 

In this approach, a family of services is described by 

— a set of states S; 

— an effect function eff : A4 x S ^ S: 

— a yield function yld : M x S* ^ {T, F, M, B}; 

— an action function act : A4 x S ^ Aau; 

satisfying the following conditions: 
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Vm€M'{3seS' yld{m, s) = M =^ Vs' e 5 • yld{m, s') {T, F}) , 



3s € S - ym G M ' 

iyld{m, s) = B A Vs' e 5 . {yld{m, s') = B ^ eff{m, s') = s)) , 

Vm e M, s e S • {yld{Tn, s) ^ M ^ act{Tn, s) = tau) . 

The sot S contains the states in which the services may be, and the functions 
eff, yld and act give, for each method m and state s, the state, reply and action, 
respectively, that result from processing m in state s. 

We define, for each s € S, a. cumulative effect function ceff ^ : M.* S m 
terms of s and eff as follows: 

ceff,{{))=s, 

ceffsia ^ (m)) = eff{m, ceff ^{a)) . 
We define, for each s e 5, a service iJg in terms of ceff yld and act as follows: 

rf{H,){a (to)) = yld{m, ceff ^{a)) , 
af{Hs){a ^ (m>) = act{m, ceff^{a)) . 

Hs is called the service with initial state s described by S*, eff, yld and act. We 
say that {Hs \ s S} is the family of services described by S, eff , yld and a,ct. 
The conditions that are imposed on S, eff, yld and act imply that, for each 
is a service indeed. It is worth mentioning that -^Hg — -ffej5f(m,s)) 
rf{Hs){{m)) = yld{m,s), and af{Hs){{m)) = act{m,s). 

7 Services for Molecular Dynamics 

In this section, wc describe a family of services which concerns molecular dy- 
namics. The formal description given here elaborates on an informal description 
of molecular dynamics given in [1]. 

The states of molecular dynamics services resemble collections of molecules 
composed of atoms and the methods of molecular dynamics services transform 
the structure of molecules like in chemical reactions. An atom can have fields 
and each of those fields can contain an atom. An atom together with the ones it 
has links to via fields can be viewed as a submolecule, and a submolecule that 
is not contained in a larger submolecule can be viewed as a molecule. Thus, the 
collection of molecules that make up a state can be viewed as a fluid. By means 
of methods, new atoms can be created, fields can be added to and removed from 
atoms, and the contents of fields of atoms can be examined and modified. A few 
methods use a spot to put an atom in or to get an atom from. By means of 
methods, the contents of spots can be compared and modified as well. Creating 
an atom is thought of as turning an element of a given set of proto-atoms into 
an atom. If there are no proto-atoms left, then atoms can no longer be created. 

It is assumed that a finite set Spot of spots such that J^,M C Spot, a total 
order < on Spot, and a finite set Field of fields have been given. It is further 
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assumed that a countable set PAtom of proto-atoms such that L ^ PAtom and 
a bijcction patom : [1, card(PAtom)] PAtom have been given. Although the set 
of proto-atoms may be infinite, there exists at any time only a finite number 
of atoms. Each of those atoms has only a finite number of fields. Molecular 
dynamics services accept the following methods: 

— for each s G Spot, a create atom method s\; 

— for each s, s' 6 Spot, a set spot method s = s'; 

— for each s, G Spot, a clear spot method s = 0; 

— for each s, s' e Spot, an equality test method s==s'; 

— for each s € Spot, an undefinedness test m,ethod .s==0; 

— for each s G Spot and v G Field, a add field method s/v; 

— for each s G Spot and v G Field, a remove field method s\v; 

— for each s G Spot and v G Field, a has field method s\v, 

— for each s G Spot and v G Field, a set field method s.v — s'; 

— for each s G Spot and v G Field, a get field method s = s'.v; 

— for each s, s' G Spot, a generate action method %s.%s'. 

Wc write A^gen for the set of all generate action methods and A^md for the set 
of all methods of molecular dynamics services. It is assumed that A4md C M. 

The states of molecular dynamics services comprise the contents of all spots, 
the fields of the existing atoms, and the contents of those fields. The methods 
accepted by molecular dynamics services can be explained as follows: 

— s!: if an atom can be created, then the contents of spot s becomes a newly 
created atom and the reply is T; otherwise, nothing changes and the reply 
is F; 

— s = s': the contents of spot s becomes the same as the contents of spot s' 
and the reply is T; 

— s = 0: the contents of spot s becomes undefined and the reply is T; 

— s==s': if the contents of spot s equals the contents of spot s', then nothing 
changes and the reply is T; otherwise, nothing changes and the reply is F; 

— s==0: if the contents of spot s is undefined, then nothing changes and the 
reply is T; otherwise, nothing changes and the reply is F; 

— s/v: if the contents of spot s is an atom and v is not yet a field of that atom, 
then V is added (with undefined contents) to the fields of that atom and the 
reply is T; otherwise, nothing changes and the reply is F; 

— s\v: if the contents of spot s is an atom and u is a field of that atom, then v 
is removed from the fields of that atom and the reply is T; otherwise, nothing 
changes and the reply is F; 

— s\v: if the contents of spot s is an atom and w is a field of that atom, then 
nothing changes and the reply is T; otherwise, nothing changes and the reply 
is F; 

— s.v = s': if the contents of spot s is an atom and v is a field of that atom, 

then the contents of that field becomes the same as the contents of spot s' 
and the reply is T; otherwise, nothing changes and the reply is F; 
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— s = s'.v: if the contents of spot s' is an atom and u is a field of that atom, 
then the contents of spot s becomes the same as the contents of that field 
and the reply is T; otherwise, nothing changes and the reply is F; 

— %s.%s': if the contents of spots s and s' are atoms and there exist /' G 
and m' £ Ai such that the contents of spot /' equals the contents of spot 
s and the contents of spot m' equals the contents of spot s', then nothing 
changes, there is no reply, and %s.%s' is turned into f.m where / and m are 
the least /' G ^ and m' & M with respect to < that satisfy the conditions 
formulated above; otherwise, %s.%s' is rejected. 

In the explanation given above, wherever wc say that the reply is T or the 
reply is F, this is meant to imply that the method concerned is turned into tau. 
Moreover, wherever we say that the contents of a spot or field becomes the same 
as the contents of another spot or field, this is meant to imply that the former 
contents becomes undefined if the latter contents is undefined. 
Let 

SS = Spot ^ (PAtom U {_L}) , 

AS = y {A^ y (^^ (PAtom U {±}))) , 

^G'Pfi„(PAtom) i^e-Pfi„(Field) 

MDS = {{a, a)GSS X AS\ rng((7) C dom(a) U {±} A 

Va G dom(a) • rng(a(a)) C dom(a) U {-L}} , 

s e MDS . 

Then wc write A4Ds for the service with initial state s described by S' = MDS U 
{t}, where t ^ MDS, and the functions eff, yld and act defined in Tables 12, 
13 and 14.^ In these tables, the functions new : 'Pfin(PAtom) — > (PAtom U {-L}), 
gencnd : SS x Spot x Spot {T, F}, and genact : SS x Spot x Spot A are 
used. These functions arc defined as follows: 

new{A) = patom{m + 1) if m < card(PAtom) , 
new{A) — 1. if m > card(PAtom) , 

where m = max{n | patom{n) G A}; 

gencnd{a, s, s') = T iff 
a{s) ^ _L A cr(s') ^ _L A 3/ G • a{f) = a{s) A 3rn G 7V4 • a(m) = a{s') ; 

genact{a, s, s') = f.m , 

® We use the following notation for functions: [] for the empty function; [d r] for 

the function / with dom(/) = {d} such that f{d) = r; f ® g for the function h with 
dom(/i) = dom(/) Udom((/) such that for all d G dom(h), h{d) = f{d) if d dom(5r) 
and h{d) — g{d) otherwise; and f ^ D for the function g with dom(p) = dom(/) \ D 
such that for all d £ dom(g'), g{d) = f{d). 
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Table 12. Effect function for molecular dynamics services 



effis\,{a,a)) = 

(cr ® [s new{dom{a))],a ® [new(dom(a)) h- > []]) if new(dom(a)) ^ -L 

eff{s\, {a, a)) — (a, a) if new(dom(a)) = ± 

e/f(6- = s', (a, a)) = (a ® [s ^ 'T(s')],a) 
effis = 0, (o-,q)) = (a® [s ±],a) 
eff{s==s', {a,a)) = {a, a) 
eff{s==0, {cr,a)) = (fT,a) 
eif (s/w, (o-,a)) = 

{a, a ® [(7(s) 1-^ a{a(s)) ® [v h-s- ±]]) if o-(s) ± A v ^ dom(a(o-(s))) 

eff{s/v, (a, a)) = (a, a) if cr(s) = _L V v € dom(Q;(cr(s))) 
eff{s\v, {a, a)) = (cr, a © [a{s) l-^ a{a{s)) {v}]) if a{s) / ± A u G dom(a(a(.s))) 

ejf{s\v, (cr, a)) = (cr, Of) if a{s) = A-W v ^ dom{a{a{s))) 
eff{s\v, (a, a)) = (cr, a) 
effis.v = s',{a,a)) = 

(a, a ® [cr(s) 1-^ a{<T{s)) ® [v h^- a-(s')]]) if o-(s) X A w € dom(a(o-(s))) 

eff{s.v = s' , (cr, a)) = (cr, a) if cr(s) = _L V v ^ dom(Q(cr(s))) 

eff{s = s'.v, {a, a)) = {a ® [s ^ a{a(s')){v)], a) if cr(s') ^ ±Av e dom(Q(a(s'))) 

eff{s = s'.v, {a, a)) = (cr, a) if ct(s') = ± V v ^ dom(a(a-(s'))) 

ejf{%s.%s' , (cr, q)) = (cr, a) if gencnd{a, s, s') — T 

ejf(%s.%s' , (cr, q)) = 1 if gencnd{a, s, s') = F 

eif (m, (cr, a)) = t if m X^d 
e#(m. [) = 1 



where / = niin{/' e \ cr{f') = <t{s)} and m = min{m' € | a{m') = cr{s)}. 
We write for the family of services described by S, eff , yld, and act. We 
write A^init for the service from this family of which the initial state is the 
unique (cr, a) G S such that dom(a) = []. 

Let (cr, a) G S, let s G Spot, let a G dom(a), and let v G dom(Q;(a)). Then 
cr(s) is the contents of spot s if cr(s) 7^ _L, is a field of atom a, and a(a)(t') is the 
contents of field v of atom a if a(a)(u) ^ _L. The contents of spot s is undefined 
if cr(s) = ±, and the contents of field v of atom a is undefined if ci;(a)(f) = ±. 
Notice that dom(a) is taken for the set of all existing atoms. Therefore, the 
contents of each spot, i.e. each element of nig{a), must be in dom(a) if the 
contents is defined. Moreover, for each existing atom a, the contents of each of 
its fields, i.e. each element of rng(Q;(a)), must be in dom(a) if the contents is 
defined. Molecular dynamics services get into state t when refusing a request to 
process a command. The molecular dynamics service with initial state t is the 
same as the divergent service D. The function new turns proto-atoms into atoms. 
After all proto-atoms have been turned into atoms, new yields _L. This can only 
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Table 13. Yield function for molecular dynamics services 



yld 


s\, (a, a)) = T 




if neui(dom(a)) 7^ _L 


yld 


s\,{a,a)) = F 




if new(dom(a)) = -L 


yld 


s = s', (o-, a)) = 


T 




yld 


s = 0, (a, a)) = 


T 




yld 


s==s', {a,a)) = 


= T 


if cr(s) = a(s') 


yld 


s==s',{a,a)) = 


= F 


if a{s) ^ o-(s') 


yld 


s=^0,{a,a)) = 


T 


if a{s) = _L 


yld 


s==0, (o-,a)) = 


F 


if cr(s) 7^ _L 


yld 


s/v,{a,a)) = T 




if a{s) ^ _L A t; dom(a(cr(s))) 


yld 


s/v,{a,a)) = F 




if cr(s) = J-W V € dom(a((T(s))) 


yld 


s\v,{a,a)) = T 




if (j(s) / -La V e dom(a(cr(s))) 


yld 


s\v, {a, a)) = F 




if cr(s) = _L V u dom(a(o-(s))) 


yld 


s\v, {a, a)) = T 




if cr(s) ^ J- Av £ dom(a(cr(s))) 


yld 


s\v, (cT,a)) = F 




if ct(s) = J-W V ^ dom(a(cT(s))) 


yld 


s.v = s' , (a, a)) 


= T 


if cr(s) / ± A w G dom(a(iT(s))) 


yld 


s.v = s' , {a, a)) 


= F 


if (t(s) = -L V V dom(a((j(s))) 


yld 


s = s' .V, {a, a)) 


= T 


if it(s') 7^ _L a V € dom(a(cr(s'))) 


yld 


s = s' .V, {a, a)) 


= F 


if fT(s') = _L V w ^ dom(a(fT(s'))) 


yld 


%s.%s',ia,a)) 


= M 


if gencnd{a, s, s') = T 


yld 


%s.%s',{a, a)) 


= B 


if gencnd{a, s,s') = F 


yld 


rn, (a. a)) = B 




if m ^ A4md 


yld 


f) = B 







Table 14. Action function for molecular dynamics services 



act{m, {a, a)) = tau if m ^ A^gen 

ac£(%s.%s', (a, a)) = genact{a, s, s') if gencnd{a, s, s') = T 

act{%s.%s' , (a, a)) = tau if gencnd{a, s, s') = F 
act{m, 1) = tau 



happen if the number of proto-atoms is finite. The initial state of AtDinit is the 
unique state in which no proto-atoms have been turned into atoms yet. 

8 Representing Programs by Molecules 

In this section, we associate each PGLD program with a PGA program for 
constructing a representation of the PGLD program by a molecule. 

Let ui ; . . . ; Uk be a PGLD program in which the foci /i, . . . , /„ and the 
methods mi, ... , m„' occur. The idea is that: 
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— an atom is created for each of the foci /i , . . . , /„ and the methods 
mi , . . . , m„' , using the; focus or method itself as the spot into which the 
atom concerned is brought on its creation, and an atom is created for each 
of the instructions ui,. . . ,Uk; 

— each atom that corresponds to an instruction is Hnked via fields to other 
atoms as follows: 

• if the corresponding instruction is of the form f.m, + f.m or —f.m, then 
the atom is linked to the atoms that correspond to the focus and method 
concerned and the atoms that correspond to the instructions with which 
execution must proceed in the cases of a positive and a negative reply; 

• if the corresponding instruction is of the form ##Z and I < I < k, then 
the atom is linked to the atom that corresponds to the instruction with 
which execution must proceed, i.e. the l-th instruction; 

— to prevent that an exception must be made of the instruction Uk in the case 
where it is of the form f.m, +f.m or — f.m, two additional atoms are created 
that are not linked to other atoms; 

— the atom that corresponds to the first instruction is brought into spot s. 

Notice that, if an atom corresponds to an instruction of the form ##Z and not 

1 < / < fc, then the atom is not linked to another atom. 

We make the assumptions that .F and M are disjoint and that PAtom is 
infinite. Under these assumptions, atom creation always leads to a positive reply 
and no test instructions are needed in the programs for constructing representa- 
tions of PGLD programs. The first assumption is only made because it permits 
this simplification. The second assumption is primarily made because there will 
be PGLD programs for which there are no PGA programs for constructing their 
representation if PAtom is finite. 

We define the PGA programs for constructing representations of PGLD pro- 
grams by means of a function pgld2md from the set of all PGLD programs to 
the set of all PGA programs. This function is defined by 

pgld2md(wi ;...; Ufe) = /i! ;...;/„!; mi! ; ... ; m„'! ; Si! ; . . . ; Sfe+2! ; 
pi(ui) ; . . . ; pk{uk) ; 
Sk+i/stop ; Sk+2/stop ; 
s = si;\ , 

where the auxiliary functions pj from the set of all primitive instructions of 
PGLD to the set of all PGA programs are defined as follows {1 < j < k): 



Pjif.m) 


= Sj /focus 


; Sj/meth ; Sj/pos ; Sj/neg ; 










Sj .focus 


= / ; Sj.m.eth = m : s j.pos = 




Sj.neg 


= Sj+l 


Pj{+f.m) 


= Sj /focus 


; Sj/meth ; Sj/pos ; Sj/neg ; 










Sj .focus 


= / ; Sj.meth = m ; s j.pos = 




Sj.neg 


= Sj+2 


Pji-f.m) 


= Sj /focus 


; Sj/meth ; Sj/pos ; Sj/neg ; 










Sj. focus 


= / ; Sj.meth = m ; Sj.pos = 


«j+2 


Sj.neg 


= Sj+l 
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Pj{##^) = Sj/jmp ;sj.jmp = Si if 1 < / < fc , 
Pj{##l) = Sj/stop if ^ (1 < Z < fc) 



and 



/i) • • • ) /n € ^ are the different foci that occur in ui ; . . . ; m^; 
mi , . . . , m„' S M are the different methods that occur in ui ; 
s, si, . . . , Sfe+2 G Spot \ U M). 



■ ■■,Uk; 



In this definition, the omitted focus is considered to be md. 

The properties stated in the following proposition show that it is easy to re- 
trieve the PGLD program P from its representation constructed by pgld2md(P) . 

Proposition 1. Let P ■ 



that MD 



(o-,a) 



ui ; . . . ; Ufc be a PGLD program, and let (a, a) he such 
pgld2md(P)| ■MDimt- Then yld{s==si, (c, a)) = T and for 



all j, I G [l,k], f G and m € M.: 

— Uj = f.m iff 

yld{u = Sj.focus, {a, a)) = T, yld{u=- 
yld{v = Sj.meth, (cr, a)) = T, yld{v=- 

yld{s = Sj.pos, (cr, a)) = T, yld{s=- 
yld{s = Sj.neg, (cr, a)) = T, yld{s=- 



+f.m iff 



=/, eff{u = Sj.focus, {a, a))) = T, 
=m, eff{u = Sj.meth, {a, a))) = T, 
=Sj+i, eff{s = Sj.pos, (cr, a))) = T, 
=Sj+i, eff{s = Sj.neg, {a, a))) = T; 

=/, eff{u = Sj.focus, (cr, a))) = T, 
=m, eff{u — Sj.meth, {a, a))) = T, 
=Sj+i, eff{s = Sj.pos, {a, a))) = T, 



yld{u ~ Sj.focus, (cr, a)) = T, yld{u= 
yld{v = Sj.meth, [a, a)) = T, yld{v- 
yld{s = Sj.pos, (cr, a)) = T, yld{s= 
yld{s = Sj.neg, {a,a)) = T, yld{s==Sj+2, eff{s = Sj.neg, {a,a))) = T; 
Uj = — f.m iff 

yld{u = Sj.focus, (cr, a)) — T, yld(u= 
yld{v = Sj.m,eth. {a. a)) = T, yld{v- 
yld{s = Sj.pos, (a, a)) — T, yld{s= 
yld{s = Sj.neg, (cr, a)) = T, yld{s= 



=/, effiu = Sj.focus, (cr, a))) = T, 
=m, eff{u = Sj.meth. (cr. a))) = T, 
=Sj+2, eff{s = Sj.pos, (a, a))) = T, 
--Sj+i, eff{s = Sj.neg, {a, a})) = T; 



— Ui 



##/ iff 



yld(s = Sj.jmp, {a, a)) = T, yld{s==si, eff{s = Sj.jmp, (cr, a))) = T; 



— Uj = for some I' ^ [l,fc] iff yld{sj\stop, (ct, a)) = T. 

Proof. From the assumption that and M are disjoint, the assumption that 
PAtom is infinite, and the definition of pgld2ind, it follows easily that: 

— all atom creations arc successful; 

— the content of each spot used in an atom creation does not change after the 
atom creation concerned; 

— Pj{uj) modifies the atom that is the contents of spot Sj only (1 < j < k), 
Sk+i/stop modifies the atom that is the contents of spot Sk+i only, and 
Sk+2/stop modifies the atom that is the contents of spot Sfe+2 only. 

Prom this, the properties stated in the proposition follow straightforwardly by 
case distinction as in the definition of Pj. □ 
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9 Interpreting Programs Represented by Molecules 

In this section, we introduce a PGA program for interpreting PGLD programs 
represented by molecules. 
The idea is that: 

— if the atom to be interpreted has the field stop, then the interpretation 
terminates; 

— if the atom to be interpreted has the field jmp, then the interpretation pro- 
ceeds with the atom to which it is linked via this field; 

— otherwise, first the basic instruction represented by the atoms to which it is 
linked via the fields focus and meth is executed and then the interpretation 
proceeds with the atom to which it is linked via the field pos or the field 
neg, depending upon the reply being positive or negative. 

The following is the PGA program for interpreting PGLD programs repre- 
sented by molecules: 

{+s\stop ; ! ; 
+s\jmp ; #9 ; 

u = s.focus ; V = s.meth ; +%u.%v ; #3 ; s = s.neg ; #4 ; s = s.pos ; #2 ; 
s = s.jmp)'^ , 

where u, w € Spot \ (.^ U M). Again, the omitted focus is considered to be md. 
Below, we write I for this PGA program. 

The program I interprets PGLD programs correctly in the sense that, for 
all PGLD programs P, the intended behaviour of P under execution coincides 
with the behaviour of the interpreter under execution on interaction with the 
molecular dynamics service that holds the representation of P constructed by 
the program pgld2md(P) when abstracted from tau. This is stated rigorously by 
Theorem 1 below. The theorem is preceded by Proposition 2 which is used in 
the proof of the theorem. The proposition concerns the local use of the spots u 
and V in /. 

For convenience, we introduce some abbreviations. Let P = ui ; . . . ; Ufc be a 
PGLD program and let i £ [l,k]. Then we write A^repr(P) for |pgld2md(P)| "md 

A^init and MDrepT{i,P) for |s = Sj ; !| •md AC'repr(P)- 

Proposition 2. Let P = ?ii ; . . . ; Ufe be a PGLD program, let i € [1, k], and let 

■^rcpr(i,P) ■SWC/I that MDT:c^,,{i^p) = |m = ; W = ; !| •md MD'j:epT(i,P)- '^^^'^ 
Ttau(|/| /md XZ?ropr(i,P)) = Ttau(|/| /md A^rcpr(<,P) ) • 

Proof. By AIP, it is sufficient to prove that for all n £ N, for all i G 
7r„(rtau(|-f| /md MDiepr{i,p))) = 7rn(rtau(|/| /md MDrepri^p))) ■ This property is 
straightforward to prove by induction on n. □ 

Theorem 1. For all PGLD programs P: 

|-P|pGLD = Ttaud-'"! /md -MVrepTiP)) ■ 
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Proof. In the proof, we make use of an auxiliary function | _ , _ | from the set of all 
natural numbers and the set of all PGLD programs to the set of all closed terms 
of sort T. It gives, for each natural number i and PGLD program ui ; . . . -jUk, a. 
closed term of sort T that denotes the behaviour of ui ; . . . ; Ufe when executed 
from position i if 1 < i < k and S otherwise. This function is defined as follows: 

\i,ui ■,...;uk\ = \tpi{ui) ; .. . ; Vfe(wfc) ; ! ; ! ; (V'i(wi) ; • ; tpk{uk) ; ! ; l)""! 

\fl<i<k , 

|i, Ui ; . . . ; Ufcl = S if 1 < i < fc 

(where tpj is as in the definition of pgld2pga). It follows easily from the definition 

of |_, _| and the axioms of PGA that if 1 < i < k: 



i. 


Ml ; . 


■ ■;uk\ 


= a o |i + 1, ux ; 


• • • 1 


Uk\ 




if Ui 


= a , 




ui ; . 


. . ; Wfel 


= |i + ; . . . 


; uk\ 


<a>\i + 2,ui; . 


• • ; Mfel 


if Ui 


= +a , 


i, 


ui ; . 


• • ; wfel 


= \i + 2,ui;...: 


; uk\ 


<a> \i + l,ui ;. 


■ ■;uk\ 


if Ui 


= -a , 


i, 


ui ; . 


• • ; Wfel 


= \l,ui;...;uk\ 








if Ui 


= ##/ 



Let P = ui ; . . . ; Ufc be a PGLD program, let 
T = {\i,P\\ie[i,k]} , 

T' = {Ttau(|/| /md A12?repr(i,P)) | i G [1, fc]} , 

and let f3 -.T ^T' be the bijection defined by 

P{\i,P\) = Ttaud-fl /md A12?repr(i,P)) • 

For each closed term p' of sort T, write (3*{p') for p' with, for all p & T, all 
occurrences of p in p' replaced by (}{p). Then, using the equations concerning 
the auxiliary function |_,_| given above and Propositions 1 and 2, it is straight- 
forward to prove that there exists a set E consisting of one derivable equation 
p = p' for each p & T such that, for all equations p = p' in E: 

— the equation /3{p) = P*{p') is derivable; 

— p' GT only if p' can be rewritten to a p" ^ T using the equations in E from 
left to right. 

Because |P|pgld = |1, f I and rtau(|/| /md MDrepT(P)) = Ttau(|-f| /md ■MDrepT{l,P)), 

this means that |-P|pgld a-nd Ttau(|-/^| /md -'^^rcpr(p)) are solutions of the same 
guarded recursive specification. Because guarded recursive specifications have 
unique solutions, it follows immediately that |P|pgld = Ttaud-'^l /md A^repr(p))- 

□ 

In the program pgld2md.(Mi ; . . . ; Uk), we could have replaced s si ; ! by 
s = Si ; Si = ; . . . ; Sk+2 = ; ! because the program / does not use the spots 
si, . . . , Sk+2- However, we cannot conceive a proof of Theorem 1 in which these 
auxiliary spots arc not found. The main problem is that wc cannot find a way 
of formulating the gist of Proposition 1, which looks to be material to a proof 
of the theorem, without referring to the spots si, . . . , Sk+2- 
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10 PGLD with Indirect Jumps 



In this section, we introduce a variant of PGLD with indirect jump instructions. 
This variant is called PGLDjj . However, preceding the introduction of FGLDij , 
we give a state-based description of the very simple family of services that con- 
stitute a register file of which the registers can contain natural numbers up to 
some bound. This register file will be used later on to describe the behaviour of 
FGLDij programs. 

It is assumed that a fixed but arbitrary number I has been given, which is 
considered the number of registers available. It is also assumed that a fixed but 
arbitrary number A'' has been given, which is considered the greatest natural 
number that can be contained in a register. 

The register file services accept the following methods: 

— for each i e [0, /] and n & [0, N], a register set method seV.r.n; 

— for each i G [0,7] and n e [0,-/V], a register test method eq:i:n. 

We write A^regs for the set {set:i:n, eq:i:n \ i G [0,1] An G [0, N]}. It is assumed 

that A^regs C M. 

The methods accepted by register file services can be explained as follows: 

— set:i:n: the contents of register i becomes n and the reply is T; 

— eq:i:n: if the contents of register i equals n, then nothing changes and the 
reply is T; otherwise nothing changes and the reply is F. 

Let s : [1,1] — > [0, A^]. Then we write Regs^ for the service with initial state 
s described hy S = {[1,1] [0,N]) U {t}, where T ^ [1,1] ^ [0,iV], and the 
functions eff, yld and act defined as follows {i e [0,1], n & [0, A''], p : [1,1] — » 
[0,7V]): 

eff{set:i:n, p) = p ® [i n] , 

ejf{eq:i:n,p) = p , 
eff{m, p) T 
effim, T) - t • 

yld{set:i:n, p) — T , 
yld{eq:i:n, p) = T 
yld{eq:i:n, p) = F 
yld{m,p) = B 
yld{m, T) = B , 

act{m,p) = tau , 
act{m, t) = tau . 

We write Regs^^n for i?e9S[i„o]e...e[/H-o]- 

We pass on to PGLD with indirect jump instructions. In PGLDy, it is as- 
sumed that there is a fixed but arbitrary finite set of foci T with regs G T and a 



if TO ^ Alregs , 



if p{i) = n , 
if p{i) i- n , 

if TO ^ Alregs , 



22 



fixed but arbitrary finite set of methods M. Moreover, we adopt the assumptions 
made about register file services above. The set {f.m \ f £ {regs},m € A4} 
is taken as the set 2t of basic instructions. 

PGLDij has the following primitive instructions: 

— for each o e 21, a plain basic instruction a; 

— for each a G 21, a positive test instruction +a; 

— for each a € 21, a negative test instruction —a; 

— for each Z e N, a direct absolute jump instruction 

— for each i G [1,/] and n G [1, A''], a register set instruction set:i:n; 

— for each i G [1,/], an indirect absolute jump instruction \i^=f^i. 

PGLDij programs have the form ui ; . . . ; Uk, where ui,. .. ,Uk are primitive 
instructions of PGLDij. 

The plain basic instructions, the positive test instructions, the negative test 
instructions, and the direct absolute jump instructions arc as in PGLD. The ef- 
fect of a register set instruction set:i:n is that the content of register i becomes n. 
The effect of an indirect absolute jump instruction is that execution pro- 

ceeds with the l-th instruction of the program concerned, where I is the content 
of register i. If is itself the l-th. instruction, then deadlock occurs. If / equals 
or Hs greater than the length of the program, termination occurs. We stipulate 
that the content of each register is initially 0. 

Like before, we define the meaning of PGLDij programs by means of a func- 
tion pgldij2pgld from the set of all PGLDij programs to the set of all PGLD 
programs. This function is defined by 

Pgldij2pgld(ui ; . . . ; Mfe) = 
V(ui);...;^K);##0;##0; 
-Fregs.eq:!:! ; ##1 ; . . . ; -|-regs.eq:l:n ; ##n ; ##0 ; 

-Fregs.eq:/:! ; ##1 ; . . . ; -|-regs.eq:/:n ; ##n ; ##0 , 

where n = min(fc, N) and the auxiliary function -0 from the set of all primitive 
instructions of PGLDij to the set of all primitive instructions of PGLD is defined 
as follows: 





= a , 






= +a , 






= -a , 






= *m 


If l<k , 




= ##0 


\f I > k , 


ip{set:i:n) 


— regs.set:i:n , 











and for each i G [1, /]: 

l^ = k + 3+{2- min{k, N) + l)-{i-l) . 
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The idea is that each indirect absolute jump can be replaced by a direct absolute 
jump to the beginning of the instruction sequence 

+regs.eq:i:l ; ##1 ; . . . ; +regs.eq:i:n ; ##n ; ##0 , 

where i is the register concerned and n = mm{k,N). The execution of this 
instruction sequence leads to the intended jump after the content of the register 
concerned has been found by a linear search. To enforce termination of the 
program after execution of its last instruction if the last instruction is a plain 
basic instruction, a positive test instruction or a negative test instruction, ; 
##0 is appended to '4>{ui) ; . . . ; V'(wfe)- Because the length of the translated 
program is greater than k, care is taken that there are no direct absolute jumps 
to instructions with a position greater than k. Obviously, the linear search for 
the content of a register can be replaced by a binary search. 

Let P be a PGLDij program. Then pgldij2pgld(P) represents the meaning 
of P as a PGLD program. The intended behaviour of P is the behaviour of 
pgldij2pgld(P) on interaction with a register file when abstracted from tau. 
That is, the behaviour of P, written IPIpoLDy, is rtau(|pgldij2pgld(P)|pGLo /regs 
RefJSinit)- 

A slightly different variant of PGLD with indirect jump instructions is in- 
troduced in [7] under the same name. 

11 The Interpretation of PGLDy Programs 

In this section, we associate each PGLDy program with a PGA program for con- 
structing a representation of the PGLDjj program by a molecule and introduce a 
PGA program for interpreting PGLDjj programs represented by molecules. This 
amounts to enhancing the PGA programs given in Sections 8 and 9. PGLDy 
programs without indirect absolute jump instructions are represented and inter- 
preted exactly as before. 

The idea of the enhancements is that: 

— an atom is created for each of the registers; 

— each atom that corresponds to a register is handled as if it concerns a di- 
rect jump instruction to the instruction with the content of the register as 
position; 

— each atom that corresponds to an indirect jump instruction is handled as 
if it concerns a direct jump instruction to the direct jump instruction that 
takes the place of the register concerned; 

— each register set instruction is handled as if it concerns an instruction for 
overwriting the direct jump instruction that takes the place of the register 
concerned. 

We define the PGA programs for constructing representations of PGLDy 
programs by means of a function pgldi j 2md from the set of all PGLDjj programs 
to the set of all PGA programs. This function is defined by 
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pgldij2md(ui ; . . . ; Ufc) = 

/i! ;...;/„!; mi! ; ... ; m„'! ; Si! ; . . . ; 5^+2! ■,s[\;...; Sj\ ; 

; . . . ; pk{uk) ; 
Sk+i/stop ; Sk+2/stop ; 

s'l/jmp ; s'l.imp = Sfe+2 ; . . . ; s^/jmp ; s'j.jmp = Sk+2 ; 

s = si ; ! , 

where the auxihary functions pj from the set of all primitive instructions of 
PGLDij to the set of all PGA programs are defined as follows {I < j < k): 



Pjif.m) 


= Sj/ focus ; Sj/meth ; Sj/pos ; 


Sj/neg 


•> 






Sj. focus = f ; Sj.meth = m ; 


Sj.pOS 




; Sj.neg = Sj+i , 


Pji+f.m) 


= Sj /focus ; Sj/meth ; Sj/pos ; Sj/neg 








Sj. focus = f ; Sj.meth = m, ; 


Sj .pos 


= Sj+l 


; Sj.neg = Sj+2 , 


Pj{-f.m) 


= Sj/ focus ; Sj/m,eth ; Sj/pos ; 


Sj/neg 


1 






Sj. focus = f ; Sj.meth = m ; 


Sj .pos 


= Sj+2 


; Sj.neg = Sj+i , 




= Sj/imp ; Sj.jmp = si 






if 1 < Z < A; , 




= Sj/ stop 






if ^(1 < / < fc) 


Pj{set:i:l) 


= Sj/reg ; Sj/ cont ; Sj/ mi ; 










Sj.reg = s[ ; Sj.cont = si ; 


.nxt = . 




if 1 < / < /c , 


Pj(set:i:l) 


= Sj/reg ; Sj/ coni ; Sj/ nxt ; 










Sj.reg = s'^ ; Sj.cont = Sk+2 '■ 




= Sj + l 


if -.(1 < / < fc) 




= Sj/jmp ; Sj.jmp = 









and 

— /i, . . . , /n € .F are the different foci that occur in mi ; . . . ; Uk', 

— mi, ... , m„/ G M are the different methods that occur in ui ; . . . ; Ufc; 

— s,si, . . . ,Sk+2,s'i, . . . ,s'j G Spot\ (J^UTW). 

The following is the PGA program for interpreting PGLDij programs repre- 
sented by molecules: 

{+s\stop ; ! ; 
+s\jmp ; #16 ; 
+s\reg ; #9 ; 

u = s.focus ; V = s.meth ; +%u.%v ; #3 ; s = s.neg ; #9 ; s = s.pos ; #7 ; 
u = s.reg ; w = s.cont ; u.jmp = v ; s = s.nxt ; #2 ; 
s = s.jmpY , 

where u,v G Spot \ (.F U A^). Below, we write I' for this PGA program. 

Theorem 2 below states rigorously that program /' interprets PGLDjj pro- 
grams correctly. In that theorem, other than in Theorem 1, we write A^repr(p) 
for |pgldij2nid(P)| •md MDinit- 
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Theorem 2. For all PGLDij programs P: 

|-P|pGLDij = Tt3ui\I'\ /md MD^cpr{P)) ■ 

Proof. The proof follows the same line as the proof of Theorem 1. This is possible 
because on interpretation any change of the state of the register file is reflected 
by a corresponding change of its molecular representation. □ 

12 Conclusions 

In this paper, we have considered the programming of an interpreter for a pro- 
gram notation that is close to existing assembly languages using PGA with 
the primitives of molecular dynamics as basic instructions. We have given PGA 
programs for constructing representations of the programs to be interpreted by 
molecules and a PGA program for interpreting those representations and we 
have shown that the latter PGA program does the interpretation correctly. We 
have experienced that, although primarily meant for explaining programming 
language features relating to the use of dynamic data structures, the collection 
of primitives of molecular dynamics in itself is suited to the programming wants 
concerned. 

We observe that: (i) the program notation for which the presented interpreter 
has been designed belongs to the simplest program notations devised ever, (ii) it 
is hard to imagine that the programs to be interpreted can be represented by 
molecules in a way that is simpler than the way chosen for the presented in- 
terpreter, (iii) it is hard to conceive of an interpreter that is simpler than the 
presented interpreter. This means not at all that the design of the interpreter was 
simple. On the contrary, the design turned out to be disappointingly difficult. It 
happened that, owing to the quest for a simple interpreter, it was inescapable 
that the design was to a great extent a trial-and-crror matter. 

Dynamic data structures modelled using molecular dynamics can straightfor- 
wardly be implemented in programming languages ranging from PASCAL [15] 
to Ctt^ [11] through pointers or references, provided that fields are not added or 
removed dynamically. Using molecular dynamics, we need not be aware of the 
existence of the pointers used for linking data. The name molecular dynamics 
refers to the molecule metaphor used in the introduction. By that, there is no 
clue in the name itself to what it stands for. To remedy this defect, we suggest 
data linkage dynamics as an alternative name. 
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